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"Abstract 

i We study the process of energy conversion from overcritical electric field into electron-positron-photon plasma. We solve 
numerically Vlasov-Boltzmann equations for pairs and photons assuming the system to be homogeneous and anisotropic. 
^ All the 2-particle QED interactions between pairs and photons are described by collision terms. We evidence several 
epochs of this energy conversion, each of them associated to a specific physical process. Firstly pair creation occurs, 
secondly back reaction results in plasma oscillations. Thirdly photons are produced by electron-positron annihilation. 
H Finally particle interactions lead to completely equilibrated thermal electron-positron-photon plasma. 
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1. Introduction 

u 

pH . Quantum electrodynamics predicts that vacuum break- 
Q^down in a strong electric field E comparable to the crit- 
' ical value E c = ro^c 3 /eh where m e is the electron mass, 
^ 'e is its charge, c is the speed of light and H is the Planck 
constant results in non-perturbative electron-positron pair 
d [production 0, 0, Nonlinear effects in high intensity 
fields can be observed already in undercritical electric field, 
y—i "see e.g. [3, [Bj]. Considerable effort has been made over 
^ last two decades in increasing the intensity of high power 
^>0 lasers in order to explore these high field regimes. Yet, the 
Schwinger field E c is far from being reached, see e.g. for 
recent review jj. There are indications that such tech- 
• nology is limited to undercritical fields due to occurrence 
£^ bf avalanches 0, H| which deplete the external field faster 
than it can potentially grow. Some authors [9| went that 
I far as to claim that "the critical QED field strength can be 
L| never attained for a pair creating electromagnetic field" . 
• While dynamical mechanisms involving increase of ini- 
tially small electric field toward its critical value appear 



problematic because of avalanches, existence of overcriti- 
cal electric field in which pair production is blocked [nj is 
widely discussed in astrophysical context in compact stars, 
e.g. hypothetical quark stars 11, neutron stars (l2| . 
see also (l3j ] for review. Pair production in such overcriti- 
cal field may occur due to several reasons, e.g. heating [3| 
or gravitational collapse of the compact object. Assuming 
existence of such overcritical electric field in this paper we 
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revisit the issue of conversion of energy from initial elec- 
tric field into electron-positron plasma once the blocking 
is released. 

The most general framework for considering the prob- 
lem of back reaction of created matter fields on initial 
strong electric field is QED. Up to now the problem has 
been treated in QED in 1+1 dimension case for both 
scalar [l5| and fermion [l6| fields. It was shown there 
that pair creation is followed by plasma oscillations due 
to back reaction of pairs on initial electric field. The re- 
sults were compared with the solutions of the relativis- 
tic Vlasov-Boltzmann equations and shown to agree very 
well. The Vlasov type kinetic equation for description of 
e + e~ plasma creation under the action of a strong elec- 
tric field was used previously e.g. in the following works 

13, m, [13. The back reaction problem in this frame- 



work was considered by [2fJ, |2l|. The kinetic theory to 



description of the vacuum quark creation under action of 
a supercritical chromo-electromagnetic field was applied 
by 22|, |23[ . Kinetic equations for electron-positron-photon 
plasma in strong electric field were obtained in [24| from 
the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy. 

Much simpler model was developed later, starting from 
Vlasov-Boltzmann equations 25| and assuming that all 



particles are in the same momentum state at a given time, 
that allowed to consider pair-photons interactions. In this 
way the system of partial integro-differential equations was 
reduced to the system of ordinary differential equations 
which was i nteg rated numerically. This model was studied 
in details in 2g,|22|, where existence of plasma oscillations 



was confirmed and extended to undercritical electric fields. 
It was also shown that photons are generated and reach 
equipartition with pairs on a time scale much longer than 
the oscillation period. 
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In this work for the first time we study the entire dy- 
namics of energy conversion from initial strong electric 
field, ending up with thermalized electron-positron-photon 
plasma which is assumed to be optically thick. With this 
goal we generalize previous treatments |15|-|27l|. In partic- 
ular, we relax the delta-function approximation of particle 
momenta adopted in j25[. In contrast, we obtain the sys- 
tem of partial integro-differential equations which is solved 
numerically on large timescales, exceeding many orders of 
magnitude several characteristic timescales of the problem 
under consideration. 

We adopt a kinetic approach in which collisions can 
be naturally described, assuming invariance under rota- 
tions around the direction of the electric field. In this 
perspective we solve numerically the relativistic Vlasov- 
Boltzmann equations with collision integrals computed 
from the exact QED matrix elements for the two parti- 
cle interactions [28| , namely electron-positron annihilation 
into two photons and its inverse process, Bhabha, Moller 
and Compton scatterings. 

The paper is organized as follows. In the next section we 
introduce the general framework, then we report our re- 
sults. Conclusions follow. Details about the computation 
are presented in the Appendix. 

2. Framework 

Based on the symmetry of the problem we consider ax- 
ially symmetric momentum space. Hence, the momen- 
tum of the particle is described by two components, one 
parallel (pii) and one orthogonal (p±) to the direction 
of the initial electric field. Then the azimuthal angle 
((/)) describes the rotation of p± around p\\. These mo- 
mentum space coordinates are defined in the following 
intervals p^ G (— oo, +00), p± G [0, +00), <fi G [0,27r]. 
Within the chosen phase space configuration, the prescrip- 
tion for the integral over the entire momentum space is 
J d 3 p — ¥ J* * d(f> dp\\ J +O ° dp± p± and the rela- 
tivistic energy is given by the following equation 



(1) 



where m is the mass of the considered particle. In the 
previous equation and hereafter we set c = h = 1 . 

In the adopted kinetic scheme the Distribution Function 
(DF) is the basic object we are dealing with and all the 
physical quantities can be extracted from it. Denoting 
with v the kind of particle, the DF /„ is commonly used 
in textbooks such that the corresponding number density 
is given by its integral over the entire momentum space 



= / d 3 p f v = 2tt 



+ OO 



d P\\ 



dp± p± fu 



(2) 



Due to the assumed axial symmetry, / does not depend on 
4> and consequently it is a function of the two components 
of the momentum only f v — fu(p\\,P±) ■ In this paper we 



use F u — 27re/„ which allows to write down Boltzmann- 
Vlasov equations in conservative form [2^, l2a |. essential 
for numerical computations. The energy density for each 
type of particle is given by its integral over the parallel 
and orthogonal component of the momentum 



dp 



+00 



dp± F v . 



(3) 



In isotropic momentum space this DF is reduced to the 
spectral energy density dp v /de. 

The time evolution of electron and positron DFs is de- 
scribed by the relativistic Boltzmann-Vlasov equation 

dF±{p\\,Pi.) , „dF±(P\\,P±) 
■ ± e hi ■ 



at 



dp\\ 



(v±(p\\,P±) ~ X±(p\\,P±) F±(P\\,P±) 



+ S(p\\,p±,E) , 



(4) 



where rf± , x± are the emission and absorption coefficients 
due to the interaction denoted by q, and the source term 
S is the rate of pair production. The sum over q covers all 
the 2-particle QED interactions considered in this work. 
In particular the electron-positron DFs in Eq. (Q}, varies 
due to the acceleration by the electric field, the creation 
of pairs due to vacuum breakdown and the interactions, 
see I Appendix A.l for details. Indeed, the Vlasov term 
describes the mean field produced by all particles, plus 
the external field. In our approach particle collisions, in- 
cluding Coulomb ones, are taken into account by collision 
terms. Particle motion between collisions is assumed to 
be subject to external field only, and the mean field is 
neglected. This is an assumption, but in dense collision 
dominated plasma such as the one considered in this paper 
this assumption is justified, see e.g. (30|. The rate of pair 
production already distributes particles in the momentum 
space according to [13| 



x log 



\eE\ 



m^(27r 



1 — exp — 



7r(mg +p 2 ± ) 
kE\ 



*(P||). (5) 



For E < E c this rate is exponentially suppressed. Besides, 
Eq. ([5]) already indicates that pairs are produced with 
orthogonal momentum, up to about m e (E/E c ) but at rest 
along the direction of the electric field. 
The Boltzmann equation for photons is 

9F ^ P±) = £ ( V ; q (P\l,Px) - X ^{p\ b Px)F 7 (p hPx ) 
q 

(6) 

and their DF changes due to the collisions only. In more 
detail, photons must be produced first by annihilating 
pairs, then they affect the electron-positron DF through 
Compton scattering. Besides, also photons annihilation 



into electron-positron pairs becomes significant at later 
times. Eqs. © and (U]) are coupled by means of the 
collision integrals, therefore they are a system of partial 
integro-differential equations that must be solved numeri- 
cally. Efficient method for solving such equations in opti- 
cally thick case was developed in |29| and later generalized 
in [28(, see |Appendix A^2] and [Appendix A. 3 for details. 

It is well known 25, 26] that both acceleration and pair 
creation terms in Eq. (jj) operate on a much shorter time- 
scale than interactions with photons described by collision 
terms in Eqs. (jl]) and ©• For this reasons we run two 
different classes of simulations, one neglecting collision in- 
tegrals which is referred to as " collisionless" and another 
one including them called " interacting" . 

3. Results 

In this section we describe our results for the collisionless 
and interacting systems separately. The boundary condi- 
tion is set when the initial electric field Eq and the initial 
DF F v q(p»,p±) are specified. For simplicity we performed 
several runs with different initial electric fields, but always 
with no particles at the beginning. When in addition to 
external electric field also particles are present from the be- 
ginning, oscillations still occur, but with higher frequency, 
as given by the plasma frequency 27J . So we start compu- 
tations with DFs null identically in the whole momentum 
space and our initial conditions can be written as 



Eq — £ E c , 
Fvo(p\\,P±) = 



£= {1,3,10,30,100} 

p± € [0, +oo) , p|| e (— oo,+oo) 



Consequently electrons and positrons are produced exclu- 
sively by the Schwinger process. 

In order to interpret meaningfully our results, we intro- 
duce first some useful quantities. Initially the energy is 
stored in the electric field and it fixes the energy budget 
available as given by 



Po 



El 



(7) 



We expect therefore the final state of the equilibrated ther- 
mal electron-positron-photon plasma to be characterized 
by the temperature 



T 




MeV , 



(8) 



where a is the Stefan-Boltzmann constant. The total en- 
ergy density of pairs p± and photons p 7 are related to the 
actual and initial electric fields, E and Eq, by the energy 
conservation law 



P±= P+ + P- 



E 2 



8 77 



Pi ■ 



(9) 



Following 2JJ, we define the maximum achievable pairs 
number density 

E 2 

Umax — 75 i (10) 



which corresponds to the case of conversion of the whole 
initial energy density into electron-positron rest energy 
density 



P±rest = (n_ +n + )m e 



(11) 



where n_ and n+ are the electrons and positrons number 
densities respectively. From the electrons and positrons 
DFs we can extrapolate their bulk parallel momentum (p||) 
as defined in Eq. (|A.5[) and the symmetry of our problem 
implies that (p\\~) = — (p\\+)- We make use of this identity 
to define the kinetic energy density of pairs 



P±ki 



P±rest 




(12) 



Therefore p±kin is the energy density as if all particles are 
put together in the momentum state with = (p\\) and 
p± = while their rest energy density is p± res t- The dif- 
ference between the total energy density and all the others 
defined above is denoted as internal energy density 



P±i 



P± ~ P±rest — P±ki 



(13) 



The term "internal" refers here to the dispersion of the 
DF around a given point with coordinates {(p\\) , (p_i_)) in 
the momentum space. 

3.1. Collisionless systems 

Since interactions with photons operate on much larger 
time scale than the pair creation by vacuum breakdown 
we first present the results obtained solving the relativistic 
Boltzmann equation Q for electrons and positrons with 
X± — V± = 0- With these assumptions we expect the 
following results to be closely related to those reported in 
27|. 




tft £ 



Figure 1: Evolution of electric field E and pairs bulk parallel momen- 
tum (p||)± obtained from the numerical solution of Eq. J4]l setting 
E = 30 E c . 




t/t c 

Figure 2: Evolution with time of the pairs energies as denned by Eqs. 
10, llTTJ, G3 and for the collisionless case E = 30 E c . All of 
them are normalized by the total initial energy po given by Eq. itTt , 

For all the explored initial conditions, there are impor- 
tant analogies between the approach adopted in [271 ] and 
the one presented in this work. For each initial field the 
first half period of the oscillation t\ is nearly equal to the 
corresponding one obtained in [27| . Also the evolution 
with time of (p\\)± during this time lapse is very similar 
to the result given by their analytic method. The time 
evolution of electric field E and (p\\)± in Compton units 
with t c — l/m e are shown in Fig. [T]for Eo = 30 E c . 

In addition to these similarities some new important fea- 
tures emerge from the current study. The manifestation of 
these new aspects is represented in Fig. [5] where we show 
how the various forms of energy defined in the previ- 
ous paragraph evolve with time. These energy densities 
are normalized to the total initial energy density po de- 
fined by Eq. ([?])• One of the most important evidences of 
this figure is that the rest energy density of pairs p±kin sat- 
urates to a small fraction of the maximum achievable one. 
This is in contrast with the result presented in [23] where 
the value given in Eq. (p~0|) was reached asymptotically. 

As a consequence the energy is mainly converted into 
other forms, namely the kinetic p±kin and internal p±i n t 
ones. Both these quantities oscillate with the same fre- 
quency but with shifted phase. Relative maxima and min- 
ima of p±kin correspond to the peaks of the bulk parallel 
momentum shown in Fig. [TJ as can be grasped from its 
definition in Eq. (fT2)l . Looking at Fig. [2]we see their rela- 
tive importance changing progressively with time. Even if 
they oscillate, the internal component dominates over the 
kinetic one as time advances. This trend points out that 
all the initial energy will be converted mostly into inter- 
nal energy, while the contribution of the kinetic one will 
eventually be small. 

In this respect, from the electron and positron DFs we 
obtain the mean squared values of the parallel and orthog- 
onal momentum as defined by Eqs. (|A.6|) and (|A.7|) . These 



Table 1: Square root of the mean squared value of orthogonal {p\)± 
and parallel (p|}± momentum, parallel momentum pp, in units of 
m e , and number density n± of pairs at the first zero of the electric 
field, saturation number density n s normalized by the maximum 
achievable one given by Eq. HOI for different initial electric fields. 



E 
E~c 




v/(pf) ± 


(Pll)i 


n x 


n s 


1 


0.4 


75 


160 


0.006 


0.018 


3 


0.8 


37 


82 


0.018 


0.037 


10 


1.3 


35 


77 


0.013 


0.041 


30 


2.0 


87 


192 


0.005 


0.016 


100 


3.5 


127 


284 


0.003 


0.011 



quantities give us some insight about the spreading of the 
DF along the parallel and orthogonal components of the 
momentum. In Tab. [T]we report their values at the end of 
runs with different initial fields. It is clear that the larger 
the initial electric field the larger is (pj_)±. This is a direct 
consequence of the rate of pair production given by Eq. 
(0 that already distributes particle along the orthogonal 
direction in the momentum space. 

The mean squared value of the parallel momentum 
(p|)± reaches a minimum value between 3 and 10 criti- 
cal electric fields. This minimum was first found in [26], 
see Fig. 3 in that paper. In Tab. [TJwe report also (p\\)±i 
which is the peak value of the bulk parallel momentum at 
the moment when electric field vanishes for the first time. 
We see from the table that also this quantity has a min- 
imum in the same range of initial fields as (pn)±- Both 
these minima are linked to the combined effects of pairs 
creation and acceleration processes. 

However, it is important to compare (pf)± and (p±)± for 
different initial fields. Indeed, this juxtaposition gives us 
quantitative informations about the anisotropy of the DFs 
in the phase space. Looking at the numerical values we 
observe how this anisotropy decreases with the increase of 
the initial electric field, which points out how an eventual 
approach toward isotropy, and therefore thermalization, 
would be much more difficult for lower initial fields. 

In Tab. [T] we compare also two different number den- 
sities n\ and n s normalized to the maximum achievable 
one. The first is the number density of pairs at the first 
zero of the electric field t%. The second is the saturation 
number density of pairs at the end of the run. We found 
that the values of ri\ are very close to the same densities 
computed in [27] with a significant amount of pairs pro- 
duced already in a very small time lapse. Let us note that 
there are maxima of both n\ and n s in the range between 
1 and 10 E c in correspondence with minima of (p|)± and 

(Pi>±- 



4 



3.2. Interacting systems 

Now we turn to the dynamics of our system on much 
larger time scales. As discussed above, in long run inter- 
actions between created pairs become important. We con- 
sider 2-particle interactions listed in Tab. IA.2l in Appendix 
and describe them by the collision integrals in Eqs. (U]) and 
([5]) using the same range of initial fields used for the colli- 
sionless systems. More sparse computational grid is used 
as calculation of collision terms imply performing multi- 
dimensional integrals in the phase space, see Appendix 
|Appendix A.2[ 

The larger the electric field the higher the rate of pairs 
production and consequently their number density. Since 
the interaction rate is proportional to particle number den- 
sities, we expect them to be important sooner for higher 
initial field. In this respect, it is worth mentioning that in 



27[ the time i 7 was estimated at which the optical depth 



for electron-positron annihilation equals unity T(i 7 ) = 1. 
There, it was found that i 7 decreases when the initial elec- 
tric field increases. Besides, the order of magnitude of 
their estimations is in agreement with the time at which 
the photons number density is around a few percent of the 
pairs number density. 

In the previous subsection dedicated to collisionless sys- 
tems, we described the anisotropy of the pairs DF by 
means of the mean squared value of parallel and orthog- 
onal momentum reported in Tab. [1] In that case, we 
knew approximately the range of orthogonal momentum 
in which the most part of electrons and positron were lo- 
cated and the orthogonal grid was chosen and kept fixed 
from the beginning. This choice was possible because the 
dispersion along the orthogonal direction was determined 
uniquely by the rate S. The extension of the grid was cho- 
sen in such a way that the value of each DF at the grid 
boundaries was small compared to the maximum value. 
For the reason that interactions redistribute particles in 
the phase space and tend to isotropize their distributions, 
the orthogonal grid must be extended to values compara- 
ble to the kinetic equilibrium temperature. To do that, we 
use initially an orthogonal grid with the same extension 
as in the collisionless system. We extend it later when 
particles are scattered toward higher orthogonal momenta 
and therefore the tails of the DF at the boundaries is not 
negligible. The extension of the parallel grid remains es- 
sentially the same as the collisionless case. 

In order to correctly describe the pairs acceleration pro- 
cess, the time step of the computation must be a small 
fraction of their oscillation period. This constraint pre- 
vents us to study the evolution up to the kinetic equilib- 
rium within a reasonable time. After hundreds of oscilla- 
tions, the energy density carried by the electric field is a 
small fraction of the pairs and photons energy densities. 
In other words, most energy has already been converted 
into electron-positron plasma. Due to this fact, the ac- 
celeration of electrons and positrons does not affect their 
DFs appreciably. This allows us to neglect the presence 
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Figure 3: Energy densities of pairs and photons obtained from the 
numerical solution of Eqs. ||4} and $§fy with initial field Eg = 100 E c . 



of the electric field hereafter. To do that we use the dis- 
tribution function at this instant as initial condition for a 
new computation in which the condition E — is imposed. 
By neglecting oscillations induced by the electric field the 
constraint on the time step of the numerical calculation 
is released, and it is now determined by the rate of the 
interactions. 

In Fig. [3] we show the time evolution of the pairs and 
photons energy densities for Eq = 100 E c . From this plot 
we can understand the hierarchy of time scales associated 
to the distinct physical phenomena we are dealing with. In 
presence of an overcritical electric field, electron-positron 
pairs start to be produced in a shortest time according 
to Eq. ([5]). As soon as they are created, electrons and 
positrons are accelerated toward opposite directions as the 
back reaction effect on the external field. The character- 
istic duration of this back reaction corresponds approxi- 
mately to the first half oscillation period. At early times, 
even after many oscillations, the energy density of photons 
is negligible compared to that of pairs, meaning that in- 
teractions do not play any role. Such a starting period, 
during which the real system can be considered truly col- 
lisionless, exists independently on the initial electric field 
even if its duration depends on it. From Fig. [3] it is clear 
that the photons energy density increases with time as a 
power law approaching the pairs energy density. 

Only when hundreds oscillations have taken place, inter- 
actions start to affect the evolution of the system appre- 
ciably and can not be neglected any further. The slope of 
the photons curve in [3] reduces indicating that pairs anni- 
hilation has become less efficient than the photons annihi- 
lation process. Now the evolution of the system is mostly 
governed by interactions. Moller, Bhabha and Compton 
scatterings give rise to momentum and energy exchange 
between electron, positron and photon populations. Be- 
sides the same collisions have the tendency to distribute 
particles more isotropically in the momentum space. Af- 
ter some time, the photons energy density becomes equal 
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Figure 4: Phase space distributions of electrons (left column) and 
photons (right column) for the initial condition Eo = 100 E c . Top: 
2.3 ■ 10 2 t c , middle: 2.3 ■ 10 2 t c , bottom: 4.6 ■ 10 6 tc- 



and then overcomes the pairs energy density. This growth 
continues until the equilibrium between pairs annihilation 

and creation processes is established e~e+ <-> 77. For this 4. Conclusions 



tial energy has been converted into rest mass energy of 
electrons and positrons. For this reason e~e + are easily 
accelerated up to relativistic velocities explaining why the 
electron DF is shifted on the right side of the phase space 
plane characterized by pn > 0. At this instant electrons are 
characterized by a relativistic bulk velocity corresponding 
to a Lorentz gamma factor 170. On the second line the 
time is 2.3 • 10 2 t c and the DFs are still anisotropic if we 
look at the parameter R introduced above. However the 
situation is different with respect to the previous stage 
because the electric field is only slightly overcritical and 
many more pairs and photons have been generated. As a 
consequence interaction rates are much larger than before 
and an efficient momentum exchange between electron and 
positron populations occurs. Both small electric field and 
collisions prevent particles to reach ultra-relativistic veloc- 
ities and for this reason the electron DF is now perfectly 
symmetric with respect to the plane p|| = 0. Only later 
on, at 4.6 • 10 6 i c for the bottom line, collisions dominate 
the evolution of the system whereas the presence of the 
electric field can be safely neglected. The pictures show 
a prominent DFs widening toward higher orthogonal mo- 
menta which is confirmed by the value R ~ 0.23. This 
remarkable evidence allows us to predict the forthcom- 
ing fate of the system to be an electron-positron-photon 
plasma in thermal equilibrium. The DFs isotropy in the 
momentum space not only indicates that the kinetic equi- 
librium condition is approached but also that the system 
is going to lose information about the initial preferential 
direction of the electric field. In the case of isotropic DF, 
the timescale on which thermal equilibrium is achieved can 
be estimated as Tth — 1/(tujt) [28|. For our anisotropic 
DF the thermalization timescale is remarkably longer. 



reason both pairs and photons curve are flat on the right of 
Fig. [3l see also [3l|. However, at this point the DF is not 
yet isotropic in the momentum space indicating that the 
kinetic equilibrium condition is not yet satisfied. In fact, 
kinetic equilibrium is achieved only at later times when 
also Moller, Bhabha and Compton scatterings are in de- 
tailed balance. At that time the electron-positron-photon 
plasma can be identified by a common temperature and 
nonzero chemical potential. This is also the last evolu- 
tion stage attainable by our study because only 2-particle 
interactions are taken into account while thermalization 
is expected to occur soon after kinetic equilibrium if also 
3-particle interactions would be included pjj ]. 

As an example, in Fig. [4] density plots of /_ and / 7 
are shown on the left and right columns respectively, for 
the initial condition Eq — 100. Their time evolution 
starts from the top line to the bottom one correspond- 
ing to three different times. After 2.3 t c both DFs are 
highly anisotropic as it is well established by the ratio 

At this stage, the electric 



R 



= 0.06. 



= ^<pjj>±/<PjjM 
field is highly overcritical and a very small fraction of ini- 



From the conceptual point of view the discovery of 
plasma oscillations [l5| as result of back reaction has been 
an important step. The next step has been the analysis of 
creation of photons from pairs |25j . In this paper we have 
studied all these phenomena in great details adopting a 
kinetic approach with two dimensional phase space. We 
have found anisotropy in momentum distributions of pairs 
and photons, which we consider one of the main results of 
this paper. Another important result is the importance of 
internal energy which limits heavily the efficiency of en- 
ergy conversion into electron-positron pairs. Both these 
effects could be in principle considered as phenomenolog- 
ical evidences of an overcritical electric field, even if they 
manifest themselves on a very short time scale. 

For the first time we studied the entire dynamics of en- 
ergy conversion from initial overcritical electric field, end- 
ing up with thermalized electron-positron-photon plasma. 
Such conversion occurs in a complicated sequence of pro- 
cesses starting with Schwinger pair production which is fol- 
lowed by oscillations of created pairs due to back-reaction 
on initial electric field, then production of photons due 
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to annihilation of pairs and finally isotropization of cre- 
ated electron-positron-photon plasma. We solved numeri- 
cally the relativistic Vlasov-Boltzmann equations for elec- 
trons, positrons and photons, with collision integrals for 
2-particle interactions computed from exact QED matrix 
elements. 

In order to appreciate the consequences of the kinetic 
treatment and the relevance of interactions separately, two 
different computations have been performed for every sin- 
gle initial condition. We called them collisionless and in- 
teracting systems in view of the fact that collision terms 
have been discarded and accounted for, respectively. The 
collisionless runs allowed us to compare our results with 
those obtained earlier, and to resolve better the momen- 
tum space of pairs. In this way we found that number 
density of pairs always saturates without exceeding 5 per 
cent of the maximum achievable number density (|10[) . in 
contrast to earlier works. This number is not far from the 
thermal number density of pairs obtained from the temper- 
ature ©. In particular, for the largest field E = 100 E c 
we obtained almost 30 per cent of the thermal number den- 
sity of pairs when the interactions are not yet important. 
It is interesting that the energy stored in initial electric 
field is mainly converted into internal and kinetic ener- 
gies of pairs, but the former becomes predominant as time 
advances. Even if the distribution in momentum space re- 
minds Maxwellian, also at the very beginning it is highly 
anisotropic, with the dispersion along the direction of elec- 
tric field exceeding orders of magnitude that in orthogonal 
direction. We conclude that simultaneous production of 
pairs and their acceleration in the same electric field is 
responsible for such peculiar form of DF of particles. 

We found that interactions become important at later 
times with respect to the average oscillation period, in 
agreement with estimates performed in [27j. For higher 
initial fields interactions become significant earlier and for 
each initial condition there is a characteristic time scale 
after which they can not be neglected. We find that pho- 
tons initially follow the distribution of pairs with nonzero 
parallel bulk momentum. 

The first equilibrium manifests itself when the perfect 
symmetry between pair annihilation and creation rates is 
established. Only later on, when scatterings have dis- 
tributed particles isotropically in the momentum space, 
the kinetic equilibrium is reached. In such state the 
electron-positron-photon plasma is generally described by 
a common temperature and nonzero chemical potentials 
for all particles and its evolution toward thermal equilib- 
rium is well understood 311 ] . 



and {1,2, ...,L — 1,L} respectively. Since we are deal- 
ing with an axially symmetric phase space with respect 
to the direction of the electric field, the parallel mo- 
mentum is aligned with it while the orthogonal compo- 
nent lays on the plane orthogonal to this preferential 
axis. Each elementary volume encloses only one momen- 
tum vector which can be written explicitly in cylindrical 
coordinates as (p\\i,p±k,<j>i)- The corresponding bound- 
aries are marked by semi-integer indices [p||i_i/2jP||i+i/2]j 
\p±k-i/2,P±k+i/2], [0i-i/2,0J+i/2]- Due to axial symme- 
try, the DFs do not depend on the azimuthal angle 4> and 
the index I will be used only to identify angles explicitly. 
We use also the symbol v which identifies the kind of par- 
ticle under consideration, {7, — , +} for photons, electrons 
and positrons respectively. From these definitions, the en- 
ergy of a particle with mass m v corresponding to the grid 
point (i, k) is 



tvik = Jrn% + pjj '. + p]_ k , (m 7 = 0, m± = m e ) . (A.l) 

In this finite difference representation the distribution 
function has a Klimontovich form and can be seen as a 
sum of Dirac deltas centered on the grid points (i, k) and 
multiplied by the energy density of particles on the same 
grid point F vik 

Fv(p\\,px) = E^H ~P\\i) S (P-L -Ptk)F vlk , (A.2) 

i k 

where Y^ik ~ Si=i Ylk=v From the definition above and 
from Eq. (J3J) the energy and number densities of particle 
v are given by 



^ ^vik 1 
ik 

~ ^ ^ T^vik 1 
ik 



(A.3) 
(A.4) 



where = F U ik/e U ik- Then the mean parallel momen- 
tum, its mean squared value and the mean squared value 
of the orthogonal momentum are 

kP\\i , (A.5) 



(A.6) 



(P\\)u = 


-E 

IK 


(P\\)u = 


ik 


(pi), = 


-E 

ik 



(A.7) 



Appendix A. Computational scheme 

The discretization of the phase space is done defining a 
finite number of elementary volumes which are uniquely 
identified by triplets of integer numbers (i, k, I). Their val- 
ues run over the ranges {1, 2, 7—1, /}, {1, 2, K—l, K} 



Due to axial symmetry the mean orthogonal momentum 
must be null identically (p±) v = . 

Appendix A.l. Acceleration and electric field evolution 

Once electrons and positrons are produced, they are ac- 
celerated by the electric field toward opposite directions. 
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The time derivative of the electron or positron parallel mo- 
mentum dpu± in the presence of an electric field E is given 
by the equation of motion 



dp 



l± 



dt 



= ±eE, 



(A.8) 



where the sign + (— ) refer to the positron (electron) and 
— e is the electron charge. 

Numerically, we move particles from one cell to another 
one such that the number of particles is conserved and Eq. 
(|A.8j) is satisfied. Acceleration causes the changing with 
time of F±ik which can written as follows 



OF. 



±ik 



dp\\ 



— ± {a>i-ik F±i-ik + ctik F±ik + oti+\k F± i+ ik ^ 

(A.9. 



where the coefficients a-^,+ are defined below 



a-i+ik 



£±l-lk P\\: 

1 



P\\i-i 



1 



P\\i~P\\i-i P\\i+i 
e±i+lfc _ P\\i 



■p\\i 



(A.10) 
(A.ll) 
(A.12) 



Also the electric field evolves according to the Maxwell 
equations. Once the currents of the moving pairs are com- 
puted, the time derivative of the electric field is known. 
Consequently a new ordinary differential equation must 
be added to the system of Eqs. dU) and ([6]). However, due 
to the uniformity and homogeneity of the physical space, 
we can describe the electric field simply using the energy 
conservation law. 

Appendix A. 2. Emission and absorption coefficients 

In kinetic theory the time derivative of the DF f v due to 
interactions between particles is generally written as 32 1 



dfu(p) 



dt 



E 

1 



(vUp)-X 9 u(p)Up)), (A.13) 



where q is the label of a specific 2-particle interaction. 
Eq. (|A.13[) represents a coupled system of partial integro- 
differential equations and can be rewritten as follows 



dF v {p) 



ot 



col! 



^J2(^ q (P)-xl(p)FAp, 



(AAA) 



where rfj 1 = 2 7r ep±_rf*. The right hand side of the previous 
equation contains the so called " collision integrals" used in 
Eqs. Q and ©. 

In order to describe how the collision integrals are com- 
puted, we write down schematically a general 2-particles 
interaction as 



1 2 

Pi P2 



3 4 

P3 P4 



which means that particles 1 and 2 having respectively 
momenta pi and P2 are absorbed; while in the same pro- 
cess particles 3 and 4 with momenta P3 and P4 are pro- 
duced. The considered interactions are shown in Table 
(|A.2p . From the kinetic theory, the absorption and emis- 



Table A. 2: Exemplification of the schematic interaction for each of 
the considered QED 2-particles interactions. 



Interaction 


1 


2 


3 


4 


Pair Annihilation 


e~ 


e+ 


7 


7 


Pair Creation 


7 


7 


e~ 


e+ 


Compton Scattering 


e ± 


7 


e ± 


7 


Bhabha Scattering 


e ± 


e=F 


e ± 




Moller Scattering 


e ± 


e ± 


e ± 





X2(P2) /2(P2) 
%(P3) 
?74(P4) 



sion coefficients for the specified process are given by 
Xi(Pi) /i(Pi) = J d 3 p 2 J d 3 p 3 J d 3 p 4 x 

x Wi i2 ;3,4 /l(Pl) /2(P2) , (A.I5) 

d 3 Pi J d 3 p 3 J d 3 Pi x 

X Wl,2;3,4 /l(Pl) /2(P2) , (A.16) 

d 3 pi J d 3 p 2 J d 3 pi x 

x 101,2:3,4 A(Pi) /a(P2) , (A.I7) 
d 3 pi J d 3 p 2 J d 3 p 3 x 

X Wl,2;3,4 /l(Pl) /2(P2) , (A.18) 



where the integrals must be calculated all over the phase 
space. The "transition rate" is given by 33] 



Wl, 2-3,4 



\M fi 



(27r) 2 16eie 2 e 3 e 4 
x 5(6! + e 2 - e 3 - e 4 ) 5 {3) (pi + p 2 p 3 - p 4 ) (A.19) 

and it contains all the informations about the probability 
that such a process occurs. The Dirac Delta's are needed 
to satisfy the energy and momentum conservation laws. 

Eqs. (|A.15[) - (|A.18I) can be discretized using the integral 
prescriptions given below 



Xliiki — 


/ 


d 3 p e p 


xi(p) /i(p) , 


(A.20) 


X2i2&2 F2i 2 k 2 


f 


d 3 p e p 


X2(p) /2(p) , 


(A.21) 


^3*3 fc 3 = 


f 


d 3 p e p 


573 (P) , 


(A.22) 


'74i 4 fe4 


f 


d 3 p e p 


Vi(p) , 


(A.23) 



where Vi n k n is the volume in the phase space which con- 
tains only the grid point with piu and p±k n - It is actually 
a ring with inner radius Pj_k-i/2> outer radius p±k+i/2 
and thickness pp+1/2 — P\\i— 1/2' An explicit expression 
for the collision integrals can be obtained inserting Eqs. 
(lXT5j) - ([X7L8)l into the corresponding Eqs. (TA30l) - (|A~23j) . 
Then we replace all the integrals over the entire momen- 
tum space with a sum of integrals over the elementary 
volumes Vik 



(A.24) 



Even if the sums are different, as well as the integrands, 
we have exactly the same sequence of integrals in all the 
emission and absorption coefficients. Due to this fact, we 
can now adopt the same treatment in order to simplify 
their expressions. 

Let us note first that the interaction cross section is in- 
variant by rotations around an arbitrary axis, therefore one 
angle can be fixed. In this respect we set <f>\ = and the 
corresponding integral gives a constant factor 2ir. Dirac 
Deltas in the transition rate wi,2;3,4 are used to eliminate 
three integrals over p^ 4 ,p±4, 4>4 and one integration over 
02 • This procedure is explained in all detail in the follow- 
ing section where the kinematics is studied. However, as 
a consequence of these choices, the momentum of the par- 
ticle 4 could differ with respect to those ones selected for 
the discrete and finite computational grid. For this reason 
Eq. (|A.20[) is no longer valid and must be modified. This 
is done " distributing" the particle 4 over three grid points 
such that number of particles, energy and momentum are 
conserved [28J. Once the correct cells are specified, we la- 
bel them with letters a, b, c and consequently the indexes 
ia,k a -,ib,kb,i c -,k c are used to identify the corresponding 
momentum components on the momentum grid. Hence, 
Eq. (IA.23P must be replaced with three equations 



V4i r k r 



/ d 3 pe p 77 4 (p), r = {a,b,c}. (A.25) 

Jv irhr 



The emission coefficients in the previous set must be mul- 
tiplied by the relative weights x a ,Xb,x c . In fact the time 
derivative of the total energy, number of particles and mo- 
mentum of the system must be null identically 



l/ik 



0. 



v i k 



t-vik 



(A.26) 
(A.27) 
(A.28) 



where the notation Q means the time derivative dQ / dt due 
to interactions only. In the previous set of equations, we 
used only the parallel momentum since the conservation 
of the orthogonal one is their direct consequence. 



The relative weights can be determined uniquely solving 
the previous system of algebraic equations. Using Dirac 
Deltas inside integrals, we can rewrite the absorption and 
emission coefficients as discrete sums 

Xlilfci -Fliifci — Eliifei ^ ' ^ ' Ri\kii 2 k 2 i 3 k 3 F2i 2 k 2 F\i\k\ > 
i 2 k 2 i 3 k 3 

(A.29) 

X2i 2 k 2 Fli 2 k 2 — ^2i 2 k 2 ^ ' ^ ' Ri\kii 2 k 2 i 3 k 3 -Fliifci F*2i 2 k 2 > 
ilki i 3 k 3 

(A.30) 

V3i 3 k 3 = e 3i 3 k 3 X! Riikii2k 2 i 3 k 3 F Ulkl F 2l2 k 2 , (A. 31) 
iifci i 2 k 2 

Vii,,k r — %r £4i r fe r ^ ^ Ri\kii 2 k 2 i 3 k 3 Fliiki ^2i 2 k 2 i 
i\k\ i 2 k 2 

(A.32) 

where the index r spans the set {a, 6, c} and the following 
coefficient is used 



R 



1 



ilkxi 2 k 2 i 3 k 3 



P-Lfca 



4 ( 27r ) 3 e u lkl e h 2 k 2 £ 3 l3 fc 3 £4 

x E J \ M M 2 ^3+1/2 - 0/3-1/2) • (A.33) 
h 

The Jacobian, due to the energy to angle change of vari- 
able, is given by Eq. ([A. 521) operating the substitution 
0i = 



J 



£i/p±k 2 



sin 4>2 (p±k 3 cos 4> h - p± kl ) - p±k 3 cos 2 sin <p h ' 

(A.34) 

Let us note that in Eqs. (IA.33[) and (|A.34I) we used €4 and 
4>2 because they are obtained from conservation laws and 
could not be associated with grid values and consequently 
labeled by indexes. Inserting Eqs. (|A.29|) - (|A.32[) into Eqs. 
(|A.26P - (IA.28|) we get a system of three algebraic equations 



^ %r E4i r fe r — £4 , 
r 

r 

^2 x rP\\i r =P\\4 



(A.35) 
(A.36) 
(A.37) 



Eqs. l|A.35p . (|A.36|) and (|A.37j) state that the global con- 
servation laws are a consequence of the conservation laws 
for each specific interaction. Now we can solve the system 

where the common denominator is 

D = P\\ a (eb - e c ) + P\\b(e c - e<t) + P\\c{e<t - e&) • (A.38) 



for the 3 unknown x r 


and its solution reads 




P\\i(eb - e c ) 4 


-P\\b(£c ~ £4) +P\\c(t4 - 


Xb = 


P\\4^c ~ e a ) H 


-PWafa ~ e c) +P\\c(£a ~ 


x c 


P\\i{ta - e b ) H 


-P\\a( e b - e 4 ) +P||b(£4 - 



Since Ri 1 k 1 i 2 k 2 i3k 3 does not depend on the DFs, it can 
be computed only once before the numerical computation 
is performed. Looking at Eqs. (|A.29I) - (|A.32|) . we easily 
see that some summations can be also computed as soon 
as the Ri 1 k 1 i 2 k 2 i3k 3 are known. Therefore we introduce the 
following quantities 



B 



^\i\k\i 2 k 2 
^2i 2 k 2 i\k\ 
Zi 3 k 3 i\k\i 2 k 2 



i 2 k 2 i 3 k 3 



or 

' Ai 3 k 3 iikii 2 k 2 



^2i 2 k 2 ^ ' Riikii 2 k 2 i 3 k 3 
i$k 3 

£313/03 Ri\k\i 2 k 2 i 3 k 3 1 
Xr ^4i r k T FLi\k\i 2 k 2 i 3 k 3 1 



and the final expression for their time derivatives are 

Xliifei -Pliifci = ^ ^\i\k\i 2 k 2 Fiiki Fi 2 k 2 i 
i 2 k 2 

X2i 2 k 2 Fli 2 k 2 = ^ ' A-2i 2 k 2 i 1 k 1 Fiiki Fi 2 k 2 i 

V3i 3 k 3 = ^ ] y ] B"Ai 3 k 3 i 1 k 1 i 2 k 2 Fi 1 k 1 Fi 2 k 2 , 
i\k\ i 2 k 2 

T \U T k T = ^ , ^ ' F>4i 3 k 3 i 1 k 1 i 2 k 2 Fi\k^Fi 2 k 2 ■ 
i\k\ i 2 k 2 

The emission coefficients requires (/ x K) 2 operations for 
each time step which must be multiplied by another factor 
(I x K) due to the computation of the analytical jacobian 
needed by the adopted method. As a result we have a total 
of about (IxK) 3 calculations at each time step which puts 
strong limits on the maximum number of grid points I and 
K. 

In order to describe processes with different timescales 
we use the Gear's method for stiff ODE's used in 28[. 
In fact, this numerical approach has an adaptive time step 
which becomes small when the DF time derivative is large, 
on the contrary it becomes large when the DF time deriva- 
tive is small. 

Appendix A. 3. Two particle kinematics 

The interaction between two particles, 1 and 2, that 
gives the particle 3 and 4 as a product can be represented 
schematically as follows 



1 + 2^3 + 4. 



(A.39) 



Each particle has 3 phase space coordinates (jp v ^ ,p v± ,4>v), 
therefore we have a total of 12 variables for the 2-particles 
interactions we are dealing with. Since each interaction 
conserves momentum and energy, they reduce to 8 inde- 
pendent degrees of freedom. That means that 4 quantities 
can be determined uniquely once the others are specified. 
For us, the 8 independent variables are px,,, px x , 0x, P2,,, 
P2 ± , P3|,, P3j_, 03; then p 4|| , p 4± , 4 , (j> 2 are functions of 
the previous ones. 



The conservation of the parallel momentum gives us the 
corresponding component of the 4-th particle 



P4„ =Pl„ +P2„ "PS,, 



(A.40) 



The orthogonal momentum for the same particle can be 
worked out using the energy conservation law 



£i + e 2 = e 3 + e 4 



(A.41) 



and using the definition of the energy given by Eq. (flj as 
follows 



(ex + e 2 - £3) 2 - Pi - m 2 



4 • 



(A.42) 



where p 4|| has been obtained in Eq. (|A.40j) . From the 
conservation of the orthogonal momentum we have the 
following relations 

PX_l COS <j)i + p 2± COS 2 = P3 ± cos 3 + p i± COS 04 (A.43) 
pi ± sin 0i + p 2± sin 2 = P3 ± sin 3 + p 4± sin 4 (A.44) 

from which we can write down analytical expressions for 
4 and 02 • Unfortunately for the previous system of equa- 
tions we have two valid solutions. For 02 we have the 
following equation 



a cos 02 + 6 sin 02 + c = , 

where the coefficients arc given by 

a = pi ± cos 0i - p 3± cos 03 , 
b = pi ± sin 0i - p 3± sin 3 , 
a 2 + b 2 + pi. - pi 



(A.45) 



(A.46) 
(A.47) 

(A.48) 



The solution of eq. (|A.45|) are given by the following con- 
ditions 

• if b ^ 0, a 2 + b 2 ^ 0, c = a 



=> 02 = -2arctan \Q 
• if b = 0, c = a =>■ 02 = 7r, 



if a ^ c, a 2 + b 2 - ac - by/a 2 + b 2 - c 2 ^ 



b'2 = 2 arctan 



b - \Ja 2 + b 2 - c 2 



if a c, a 2 



ac 



+ b\/a 2 + b 2 - c 2 



62 = 2 arctan 



b + \Ja 2 + b 2 - c 2 



a — c 



Once 02 has been chosen, 4 can be easily obtained from 
the Eqs. (|A~43|) and (|A~44|> . 

The Jacobian of Eq. (|A.34I) has been computed using 
the following identity for the Dirac Delta 



*(/(*)) = £ 



8(x - Xj) 
\(df/dx) Xi \ ' 



(A.49) 
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where / is a function such that f[Xi) = 0. In our frame- 
work the function inside the Dirac delta is given by the 
energy conservation 



f(h) = e i + £ 2 - £3 ~ £4(^2) ■ 



(A.50) 



where 4>2 is now the independent variable. From the previ- 
ous equation we compute its derivative with respect to </>2 
and the value ^2 such that /(^j-j) = 0. Rewriting explicitly 
Eq. (|A.49j) we have that 

<W 2 » = } | = J5 ^ - ^) . (A.51) 

where the explicit equation for J is 

Q 



J 



p 2± |sin(0 2 - 4>a)\ 



(A.52) 



Acknowledgments 

AB is supported by the Erasmus Mundus Joint Doc- 
torate Program by Grant Number 2010-1816 from the 
EACEA of the European Commission. 

References 

[1] F. Sauter, Uber das Verhalten eines Elektrons im homoge- 
nen elektrischcn Feld nach der relativistischen Theorie Diracs., 
Zeits. Phys. 69 (1931) 742-+. 

[2] W. Heisenberg, H. Euler, Zeits. Phys. 98 (1935) 714 — h. 

[3] J. Schwinger, On Gauge Invariance and Vacuum Polarization, 
Phys. Rev. 82 (1951) 664-679. 

[4] C. Bula, K. T. McDonald, E. J. Prebys, C. Bamber, S. Boege, 
T. Kotseroglou, A. C. Melissinos, D. D. Meyerhofer, W. Ragg, 
D. L. Burke, R. C. Field, G. Horton-Smith, A. C. Odian, J. E. 
Spencer, D. Walz, S. C. Berridge, W. M. Bugg, K. Shmakov, 
A. W. Weidemann, Observation of Nonlinear Effects in Comp- 
ton Scattering, Physical Review Letters 76 (1996) 3116—3119. 

[5] D. L. Burke, R. C. Field, G. Horton-Smith, J. E. Spencer, 
D. Walz, S. C. Berridge, W. M. Bugg, K. Shmakov, A. W. Wei- 
demann, C. Bula, K. T. McDonald, E. J. Prebys, C. Bamber, 
S. J. Boege, T. Koffas, T. Kotseroglou, A. C. Melissinos, D. D. 
Meyerhofer, D. A. Reis, W. Ragg, Positron Production in Mul- 
tiphoton Light-by-Light Scattering, Physical Review Letters 79 
(1997) 1626-1629. 

[6] A. Di Piazza, C. Miiller, K. Z. Hatsagortsyan, C. H. Keitel, 
Extremely high-intensity laser interactions with fundamental 
quantum systems, Reviews of Modern Physics 84 (2012) 1177— 
1228. larXiv: 11117 3886 doi : 10 . 1103/RevModPhys . 84 . 1177] 

[7] A. R. Bell, J. G. Kirk, Possibility of Prolific Pair 
Production with High-Power Lasers, Physical Review 
Letters 101 (20) (2008) 200403. arXiv:0808.2107 
|dol : 10 . 1103/PhysRevLett ■ 101 . 200403] 

[8] S. S. Bulanov, T. Z. Esirkepov, A. G. R. Thomas, J. K. Koga, 
S. V. Bulanov, Schwinger Limit Attainability with Extreme 
Power Lasers, Ph ysical Review Letters 105 (22) (201 0) 220407. 
larXiv:1007.4306l |dol:10. 1 103/PhysRevLett . 105 . 220407 

[9] A. M. Fedotov, N. B. Narozhny, G. Mourou, G. Korn, Limita- 
tions on the Attainable Intensity of High Power Lasers, Phys- 
ical Review Letters 105 (8) (2010) 080402. a rXiv: 1004 .5398 
|dol : lb~1 103/PhysRe vLett . 105 . 080402] 
[10] V. V. Usov, Bare Quark Matter Surfaces of Strange Stars 
and e+e~ Emission, Phys. Re v. Lett. 80 (1998) 230-233. 
arXiv:arXiv:astro-ph/9712304 



[11] C. Alcoc k, E. Farhi, A. Olinto , Strange stars, ApJ 310 (1986) 
261-272. |doi: 10. 1086/1646791 

[12] R. Belvedere, D. Pugliese, J. A. Rueda, R. Ruffini, S.-S. Xue, 
Neutron star equilibrium configurations within a fully rela- 
tivistic theory with strong, weak, electromagnetic, and grav- 
itational interactions, Nuclear Physics A 883 (2012) 1-24. 
larXiv: 1202 . 65001 [doi : 10 . 10 16/ j .nuclphysa. 2012 T02T7jl8"| 

[13] R. Ruffini, G. Vereshchagin, S.-S. Xue, Electron-positron pairs 
in physics and astrophysics, Physics Reports 487 (2010) 1—140. 

[14] V. V. Usov, Thermal Emission from Bare Quark Matter 
Surfaces of Hot Strange Stars, ApJ 550 (2001) L17 9-L182. 
arXiv:arXiv:astro-ph/0103361 , doi : 10 . 1086/319639 

[15] Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, E. Mottola, 
Pair production in a strong electric field, Phys. Rev. Lett. 67 
(1991) 2427-2430. |doi : 10 . 110 3/PhysRe vLett . 67 ■ 2427] 

[16] Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, 
E. Mottola, Fermion pair production in a strong 
electric field, Phys. Rev. D45 (1992) 4659-4671. 
|doi : 10 . 1 103/PhysRevD ■ 45 ■ 4659 1 

[17] T. Matsui, K. Kajantie, Decay of strong color electric field 
and thermalization in ultra-rclativistic nucleus-nucleus colli- 
sions, Physics Letter B 164 (1985) 373. 

[18] S. A. Smolyansky, G. Roepke, S. Schmidt, D. Blaschke, V. D. 
Toneev, A. V. Prozorkevich, Dynamical derivation of a quan- 
tum kinetic equation for particle production in the Schwinger 
mechanism, ArXiv High Energy Physics - Phenomenology e- 
printf.arXiv : arXiv : hep- ph/97 12377 

[19] S. Schmidt, D. Blaschke, G. Ropke, S. A. Smolyan- 
sky, A. V. Prozorkevich, V. D. Toneev, A Quan- 
tum Kinetic Equation for Particle Production in the 
Schwinger Mechanism, International Journal of Modern 
Phys ics E 7 (1998) 709-722. |arXiv:arXiv7i iep-ph/9809227 
doi : 10 . 1142/S0218301398000403 

[20] G. Gatoff, A. K. Kerman, T. Matsui, Flux-tube model for 
ultrarclativistic heavy-ion collisions: Electrohydrodynamics 
of a quark-gluon plasma, Phy s. Rev. D36 (1987) 114-129. 
doi : 10 . 1 103/Phy sRevD .36.114 

[21] Y. Kluger, E. Mottola, J. M. Eisenberg, Quantum Vlasov 
equation and its Markov limit, Phys. Rev. D58 (12) 

(1998) 125015. |arXiv : arXiv : hep-ph/9803372| 

doi : 10 . 1103/PhysRevD . 58 . 125015 

[22] J. C. R. Bloch, V. A. Mizerny, A. V. Prozorkevich, C. D. 
Roberts, S. M. Schmidt, S. A. Smolyansky, D. V. Vinnik, Pair 
creation: Back reactions and damping, Phys. Rev. D60 (11) 

(1999) 116011. larXiv : arXiv : nucl-th/9907027 1 

doi: 10. 1 103/Phy sRevD. 60. 116011 

[23] D. V. Vinnik, A. V. Prozorkevich, S. A. Smolyansky, 
V. D. Toneev, M. B. Hecht, C. D. Roberts, S. M. 
Schmidt, Plasma production and thermalisation in a strong 
field, European Physical Journal C 22 (2001) 341-349 . 
larXiv : arXiv : micl- th/0103073 doi : 10 . 1007/sl00520100787 

[24] D. B. Blaschke, V. V. Dmitriev, G. Ropke, S. A. Smolyansky, 
BBGKY kinetic approach for an e — e+7 plasma created from 
the vacuum in a strong laser-generated electric field: The one- 
photon annihilati on channel, Phys. Rev. D84 (8) (2011 ) 085028. 
arXiv : 1 105 . 5397 doi : 10 . 1 103/PhysRevD . 84 . 085028 

[25] R. Ruffini, L. Vitagliano, S.-S. Xue, On plasma oscillations 
in strong electric fields, Phys. Lett. B 559 (2003) 12-19. 
|arXiv:ar Xiv : astro-ph/03025l9~| 

[26] R. Ruffini, G. V. Vereshchagin, S.-S. Xue, Vacuum polarization 
and plasma oscillations, Phys. Lett. A 371 (2007) 399-405. 

[27] A. Benedetti, W.-B. Han, R. Ruffini, G. V. Vereshchagin, On 
the frequency of oscillations in the pair plasma generated by 
a strong electric field, Physics Letters B 698 (2011) 75 -79. 
arXiv : 1102 . 4287] |doi : 10 . 1016/ j .physletb . 2011 . 02 ■ 050| 

[28] A. G. Aksenov, R. Ruffini, G. V. Vereshchagin, Thermalization 
of the mildly relativistic plasma, Phys. Rev. D79 (4) (2009) 
043008. arXiv : 0901 . 4837] [doi : 10 . 1103/PhysRevD . 79"T0 43008 

[29] A. G. Aksienov, M. Milgrom, V. V. Usov, Structure of Pair 
Winds from Compact Objects with Application to Emis- 
sion from Hot Bare Strange Stars, ApJ 609 (2004) 363-377. 



11 



|arXiv : arXiv : astro-ph/0309014) |doi : 10 . 1086/421006 

[30] S. R. de Groot, W. A. van Leeuwen, C. G. van Weertn, Rel- 
ativistic kinetic theory. Principles and Applications, North- 
Holland, Amsterdam, 1980. 

[31] A. G. Aksenov, R. RufBni, G. V. Vereshchagin, Thermal- 
ization of Nonequilibrium Electron-Positron-Photon Plasmas, 
Phys. Rev. Lett. 99 (12) (2007) 125003-+. 

[32] D. Mihalas, B. W. Mihalas, Foundations of Radiation Hydro- 
dynamics, New York, Oxford University Press, 1984. 

[33] L. D. Landau, E. M. Lifshitz, Physical Kinetics, Elsevier, 1981. 



12 



